All-optical Realization of Quantum Ratchets 
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A theoretical realization of an all-optical quantum ratchet is proposed in a medium composed of an array of 
coupled waveguides. By coupling light into two adjacent waveguides, and calculating the expectation value 
for the position space operator, we demonstrate the quantum ratchet like behavior of this optical system. © 
2013 Optical Society of America 



A ratchet is a device in which there is directed motion 
of a particle in one direction, but in which motion is 
blocked in the opposite direction. The Brownian mo- 
tion ratchet has attracted considerable attention since 
the ratchet paradox was resolved by Feynman [1]. The 
classical ratchet, which has been studied extensively in 
the context of biological motors, requires a combination 
of an asymmetric potential and a dissipative, random 
noise process, which together are responsible for the di- 
rected ratchet motion [2]. In recent years, a quantum 
version of the ratchet has been demonstrated in Bose- 
Einstein condensates. Unlike the classical ratchet which 
requires a symmetric input and an asymmetric potential, 
the quantum ratchet does not require a dissipative pro- 
cess, and one can utilize a symmetric potential with an 
asymmetric input. The directed motion is a result of an 
underlying asymmetry between the potential in which 
the particle is moving and the quantum mechanical den- 
sity distribution. This is illustrated in Fig. 1. In part 
(a), there is a symmetry between the potential and the 
density distribution, and there is no net motion of the 
particles because when the potential is ON, the parti- 
cle is always trapped at the minima of the potential. In 
part (b), the asymmetry between the potential and the 
density distribution leads to net motion of the particle. 
A standard model to describe quantum ratchets is the 
kicked rotor given by the Hamiltonian, 
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where the periodic potential is produced by a standing 
laser field. This field is applied at t = nT where n is an 
integer and the field consists of a series of short pulses. 
The initial state, ip, of the particle is a combination of 
states with different momenta, 



i/j(x) = ip (e iqx + ae i(p e iqx+2ikx ) 



(2) 



The phase (p and the parameter a introduce the asym- 
metry between the potential and the density distribution 
and lead to directed motion of the particles. This has 
been extensively studied both theoretically and experi- 
mentally [4]. In what follows, we suggest an all-optical 
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Fig. 1. (Color Online) A rachet results from the asym- 
metry of the potential and initial density distribution. 
If the potential (black) and the distribution (blue) are 
symmetric (a), there is no net motion of the particles 
(gray and black circles). However, if the potential is not 
symmetrical with respect to the distribution (b), there 
is a net flow of particles to the left; following [3]. 



realization of quantum ratchets by using coupled array 
of waveguides. A previous, related work, has investigated 
the optical ratchet effect in coupled waveguide arrays [5] . 
We note that such arrays of waveguides have proved 
to be of great utility in the realization of a number 
of effects from condensed matter physics and quantum 
physics, such as Bloch oscillations [6], Anderson localiza- 
tion [7,8], quantum Zeno effect [9], electromagnetically 
induced transparency [10], Dirac zitterbewegung [11,12], 
Talbot effect [13], quantum random walks [14], surface 
solitons [15, 16], Zener tunneling [17], and a quantum 
bouncing ball [18], among others. 

Consider an array of waveguides, shown in Fig. 2, in 
which the neighboring waveguides are evanescently cou- 
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pled. This system is described by a tight binding Hamil- 
tonian of the form, 

H = h ^ p(j)a] aj + HC £>t +1 aj + ot_ l0i ) (3) 



where the array is labeled by the index j and (3(j) is 
related to the refractive index of the jth waveguide. In 
Eq. (3), C is the coupling between adjacent waveguides, 
and Oj (aj) is the creation (annihilation) operator for the 
jth waveguide. These obey the Heisenberg equations, 



-iP(j)dj - iC(a j+1 + dj-x) 



(4) 



From now onwards, we consider a classical description 
(for quantum mechanical descriptions, see [19]). Thus, 
the a/s would be regarded as numbers. 
We now write 
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a(k)e ikj dk 



dj = ■ a 

3 ' 
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We can then write the equations in the Fourier space 
using Eq.(4) and Eq.(5), as 



a(k) = —2iCa{k) cos k — i/3 



d 



d(-ik) 



a(k) (7) 



This should now be compared with the Hamiltonian in 
Eq. (1). In Eq. (7), we already have the periodic poten- 
tial. In Eq. (1), we have a term that is quadratic in p, i.e. 
a term like . However, for a Bose-Einstein conden- 
sate (BEC) [20], the distribution in p is very narrow and 
thus an expansion around some mean momentum would 
reduce p 2 to a term linear in p. On the other hand, for 
the waveguide structure we can make the choice, 
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and then Eq.(7) reduces to 



a(k) = —2iCa(k) cos k — U^^7 — ^y^(^) ^ 

The choice of the waveguide index, Eq. (8), has been 
experimentally realized in studies on Bloch oscillations 
in waveguide structures. It should be borne in mind that 
the roles of position and momentum have been reversed 
in the optical realization of quantum ratchets. In Eq. (1), 
the periodic potential is in coordinate space whereas for 
the waveguide structure the periodic potential is in the 
Fourier space, k. In the quantum case, one studies the 
mean values of p, i.e. < p >, to obtain directed motion. 
In our waveguide system, we have to study the behavior 
in the site space, j. In the quantum case, we specify the 
wavefunction, Eq. (2), in the coordinate space. Here, the 
corresponding analog would be the value of a(k) in the 
k space, which can be translated to the site space j. The 



appropriate input conditions on the field amplitude will 
be 



a,j(t = 0) = Sj^o + aSj^e 2 



(10) 



We assume that the array index is from j = — oo to oo, 
see Fig. 2, where the j = is the middle waveguide. 



The 



play the role of the density distribution. The 



analogs of quantum ratchets that we would find would 
be due to a ^ 0, ip ^ 0. The familiar Bloch oscillations 
occur for a = 0. The possibility of additional interference 
effects due toa^O leads to the realization of the analog 
of quantum optical ratchets. 
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Fig. 2. A large waveguide array with the input fields 
shown. 

In order to study the ratchet problem, we calculate 
the expectation values for < j > and < j 2 > which are 
analogous to the expectation values for momentum and 
energy respectively. These expectation values are calcu- 
lated as < j >= Y%L-oo3lj and < f >= T^-oo^h 
where Ij is the output intensity of the j th waveguide. We 
write the solution to the coupled equations, Eq. (4), in 
the form 



(ii) 



where dj(0) is given by Eq. (10) and the Green's function 
is given by Eq. (21) from [21] 
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Using the initial condition of Eq. (10), the output inten- 
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sity from the waveguide array can be written as 

Ij = \Gj,o\ 2 + \aG jA \ 2 +aG jfi Gl ie iv +aG* fi G 3A e iv 
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(13) 



Fig. 3(a) shows the output intensity from the waveg- 
uide array, given by Eq. (13), when a = 0. We note 
immediately from Fig. 3(b) that there is an asymmetry 
in the profile along the j-axis, and this asymmetry is ev- 
idence for directed motion. The case when a = 0, in Fig. 
3(a), results in symmetric profiles, and corresponds to 
the well-known Bloch oscillations [6,21]. 
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Fig. 3. (Color Online) Output intensity distribution, ij, 
as a function of z for j = — 12, 12. The values of the 
parameters are /3/C = 0.73, cp = 37°, and (a) a = 0, (b) 
a = l. 



result of the interference since the second term in both 
Eq. (14) and(15) are identical. Also, the dependence of 
energy on the phase is identical to that of the momentum 
as seen in Fig. 4(b) and 5(b). 
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Fig. 4. (Color Online) Plots showing ratchet like behavior 
for (a) the momentum and (b) the energy. The values of 
the parameters are a = 1, /3/C = 0.73, and cp = 37°. 
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Fig. 5. (Colore Online) Plots showing ratchet like behav- 
ior for (a) the momentum and (b) the energy. The values 
of the parameters are a = 1, j3/C = 0.73, and cp = 217°. 



Using Eq. (13), the expectation value for the site po- 
sition can be written as (after using properties of Bessel 
functions) 

E^ = H 2 + ^sin(|)sin(|-^. (14) 

j = — OO 

Similarly, the expectation value for the analog of energy 
is given by 

V* -2r i ,2 4aC . ( f3z\ . ( f3z \ 

j— — OO 7 

1 + M 2 (AC . (I3z\\ 2 

For small values of z, the average position, Eq. (14), is 
proportional to — asm(cp) z, which shows that direction 
of transport for the photons is dependent on the the 
relative phase. Fig. 4(a) shows that the photons with a 
phase of 37° are directed to the left in Fig. 2 for small 
values of z while Fig. 5(a) shows the photons with a 
phase of 217° are directed to the right. The directed 
motion is a result of the interference of the two photons 
as the second term in Eq. (14) is linear in a. The analog 
of the energy is also linear in z for small values of z as a 



In summary, we have described an all-optical real- 
ization of the quantum ratchet in a system of evanes- 
cent ly coupled waveguide arrays. Starting from the tight- 
binding Hamiltonian that is used to describe such a 
system, one can derive an equation of motion that is 
analagous to the kicked rotor model of the quantum 
ratchet. The principal difference between the theoreti- 
cal frameworks is that in our optical realization of the 
quantum ratchet, the periodic potential is in the Fourier 
space instead of coordinate space. This requires one to 
study the optical ratchet behavior in the site-space. This 
is fortunate because all waveguide studies are done in 
site-space. We have demonstrated that by choosing a 
suitable input to the waveguide array that consists of 
a linear combination as described by Eq. (10), one can 
recover the Bloch oscillations as a limiting case of the 
ratchet behavior when a = 0. 
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